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Abstract Purpose: We are building an MRTOF-MS 
(Multi-Reflection Time-Of-Flight Mass Spectrometer) 
for isobaric separation for the Lanzhou Penning Trap. 
The potentials applied on the electrodes of our MRTOF 
mass analyzer operating in in-trap-lift mode have to be 
optimized to achieve a very high mass resolving power. 

Methods: Our method to design and optimize an 
MRTOF mass analyzer has been updated to introduce 
constraints on the potentials and this method now can 
be used to optimize the parameters of MRTOF-MS 
both operating in mirror-switching mode and in in- 
trap-lift mode. By using this method, the optimal po- 
tential parameters of the electrodes have been obtained 
for our MRTOF mass analyzer operating in the in-trap- 
lift mode. 

Results and Conclusion: With a beam size of 2.8 mm 
diameter and an initial average ion kinetic energy of 
1500 eV, the maximal mass resolving power has been 
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achieved to be 3.2 x 104 with a total TOF of 7.0 ms for 
an ion species of 4°Ar'*. It can reach up to 5.6 x 104 
for a beam size of 0.3 mm diameter. The simulation 
shows that the inaccuracy of the potentials applied on 
the outermost mirror electrodes M1-M2 must be less 
than 50 ppm or preferably 20 ppm. 


Keywords Time-of-flight mass spectrometer - 
Multiple-reflection - Mass measurement - Isobaric 
separation - Exotic nuclei 


1 Introduction 


The Multi-Reflection Time-Of-Flight Mass Spectrom- 
eter (MRTOF-MS) has been developed as a new de- 
vice in recent years. By extending the flight path us- 
ing multi-reflection between electrostatic ion mirrors, 
an MRTOF-MS can reach a very high mass resolv- 
ing power of >100,000 in a compact structure. More- 
over, it also has other unique advantages, such as ex- 
tremely short measurement time, a large mass range, 
very high sensitivity and non-scanning operation. Plaf 
et al. [1] reviewed the MRTOF-MS developments for 
the research with short-lived nuclei and different instru- 
mental implementations. Up to now, many MRTOF- 
MSs for mass measurements and isobaric separations 
have been commissioned [2-6] or under construction [7- 
11]. Many good results have been achieved, as examples, 
the masses of 8?Zn [12], °3:54Ca [13] and °?5°K [14] have 
been measured by using MRTOF-MS at ISOLDE/CERN 
and the masses of heavy nuclides produced by fusion- 
evaporation reactions at GARIS-II/RIKEN [15]. 
Injection and ejection of ions into/from MRTOF- 
MS have been achieved by switching the electric po- 
tentials of the electrodes to appropriate values while 
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the ions are passing. Thus, according to the switch- 
ing mode, all MRTOF-MSs can be cataloged into two 
types: mirror-switching and in-trap-lift modes. In the 
former mode, the electric potentials of the entrance and 
exit mirrors are switched to lower values and the po- 
tential of the intermediate drift-tube keeps at 0, while 
in the latter, all the potentials of the mirrors keep at 
their optimal values and the intermediate drift-tube is 
switched to its appropriate value, in other words, the 
drift-tube acts as a lift to the ions and the energies of 
the ions vary, thus the ions can be injected, trapped 
and ejected by only switching the potential of the drift- 
tube. Wolf et al. [16] gave a detail theoretical study, 
and many MRTOF-MSs [4,6,10] have adopted the in- 
trap-lift mode. 

The Lanzhou Penning Trap (LPT) [17], of which the 
main task is to perform direct mass measurements on 
atoms and nuclei with high precision, is presently under 
construction at the Institute of Modern Physics, Chi- 
nese Academy of Sciences. We are building an MRTOF- 
MS for isobaric separation for LPT. The design and 
optimization procedures have been reported in detail 
in Ref. [7,18], where the mirror-switching mode has 
been assumed and a maximal resolving power has been 
achieved to be 1.3 x 10° with a total time-of-flight of 
6.5 ms for an ion species of *°Ar!*. In this paper, we 
report the optimization of our MRTOF mass analyzer 
operating in in-trap-lift mode and the corresponding 
results. 


2 Parameters of the MRTOF mass analyzer 


Fig. 1 shows the geometry of our MRTOF mass ana- 
lyzer. It consists of two identical electrostatic mirrors 
(each contains a set of four electrodes) in combination 
with two einzel lens and a drift-tube. All electrodes of 
this analyzer have a cylindrical shape. The total length 
is 708 mm with an inner diameter of 60 mm. The mirror 
electrodes M1-M4 (numbered from the outermost) have 
lengths of 20, 16, 26 and 26 mm, respectively, and the 
lens electrode L has 46 mm. The intermediate drift-tube 
D has a length of 400 mm. The distance between two ad- 
jacent electrodes is 4 mm. The time focus plane where 
places a micro-channel-plate detector or a Bradbury- 
Nielsen gate [19] is set to locate 255 mm after the ana- 
lyzer. Different from our previous optimization, we set 
the mass analyzer operating in the in-trap-lift mode. 


3 Update of the optimization procedure 


In our previous paper [7], we reported an optimization 
strategy in detail to find the optimal parameters for our 
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Fig. 1 The geometry of our MRTOF mass analyzer. 


mass analyzer operating in the mirror-switching mode. 
It includes two sub-procedures, a global search and a 
local refinement. The ion trajectories are calculated in 
the SIMION [20] according to our specified user pro- 
gram as usual, and the variations of the parameters are 
performed in a separate optimization program that is 
coded by ourselves by using the Nelder-Mead simplex 
algorithm [21]. Because we search all the possibilities in 
the whole parameter space the best parameters can be 
found for sure. 

We found the drawback of our optimization pro- 
gram when we tried to optimize our MRTOF mass an- 
alyzer operating in the in-trap-lift mode. Practically, 
the kinetic energy of the ions from the RFQ cooler and 
buncher in the upstream has to be fixed at a certain 
value, say, 1500 eV. It constrains all the potentials ap- 
plied on the electrodes of the mass analyzer at <1500 V. 
Thus, the optimization has to be done according to 
boundaries on the possible potentials on the electrodes. 
Our previous program cannot do this and we have to 
update it. 

The Nelder-Mead simplex algorithm originally does 
not support constraints (upper and lower boundaries). 
A penalty function or a domain transformation was pro- 
posed to use. Box [22] proposed a way to handle some 
explicit and some implicit constraints and Guin [23] 
later refined the idea. Le Floc’h [24] compared all these 
methods. To solve our problem, we perform the follow- 
ing strategy. When a trial point is out of the bound- 
aries and the centroid’s coordinate P; respects the con- 


straints, the point is reset at l + w(P; — l) if the lower 
boundary | is broken or at h — w(P; — 1) if the upper 
boundary h is broken, where w is a uniform random 
number between 0.00001 and 0.5. If the centroid breaks 
the constraints, the point is reset at l + 0.00001 if the 
lower boundary is broken or at h — 0.00001 if the upper 


boundary is broken. 


4 Results and discussion 


To calculate the mass resolving power with respect to 
the number of revolutions, we considered the following 
as the initial beam conditions. The beam consists of 
100 ions of 4°Ar!*, of which the average kinetic energy 
is 1500 eV, and the bunch width is 20 ns (FWHM) 
at the middle of the analyzer. (Since the calculation 
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is very time-consuming we cannot increase the num- 
ber of ions due to our existing computation capacity.) 
In the middle between the ion mirrors the ion beam 
coordinates x and y orthogonal to the optical axis z, 
the angles a=dx/dz and B=dy/dz with respect to the 
axis z and the energy K are approximately represented 
by Gaussian distributions with the standard deviations 
O,=0y=1.0 mm, og=og=1.5 mrad and og=8.5 eV. 
All these parameters are the same as those in Ref. [7]. 
For comparison, we also reduce the beam spot from 
0,=0y=1.0 mm to 0.5, and 0.1 mm, corresponding to 
beam sizes with diameters of 2.8, 1.4 and 0.3 mm, re- 
spectively. If not specified, the results are obtained for 
the beam size with 2.8 mm diameter. 

The mass resolving power R, which is one of the 
most important quantities in mass spectrometry, is cal- 
culated as R = TOF /2/ATOF, where TOF and the 
overall spread ATOF are the centroid and the full- 
width-at-half-maximum (FWHM) of the time-of-flight 
peak. 


4.1 Optimized results for the beam size with 2.8 mm 
diameter 


Considering the symmetric shape of our MRTOF mass 
analyzer and the practical potential controls through 
the power supplies during the experiments, we set the 
potentials applied on the four mirror electrodes, one 
lens electrode and the intermediate drift-tube as the 
optimization parameters, six in total. Because the in- 
coming ions have an average energy of 1500 eV, all the 
potentials applied on the electrodes have been set an 
upper limit of 1400 V. 

In the global search step, we totally considered 4 x 
4x4x7x9x 3 = 12096 potential combinations as the 
initial parameter sets and at last chose only the best 
five new sets obtained from the search for further local 
refinement. By the Nelder-Mead simplex algorithm all 
parameters varied and a large number of local minima 
were obtained. We chose the parameter set which gives 
the biggest value of mass resolving power as the op- 
timal one, rather than the smallest peak width of the 
TOF when the mass analyzer operates in the mirror- 
switching mode; because in the in-trap-lift mode the 
average kinetic energy of the ions travelling in the trap 
varies according to the formula, Kintrap = e[1500 — 
Uari ft], Where Uarift is the potential of the drift tube, 
and the total TOF varies with different parameter set. 

We optimized the potentials for different number of 
revolutions and obtained the maximal resolving power 
that can be achieved, as shown in Fig. 2. For the beam 
size with 2.8 mm diameter, when increasing the number 
of revolutions N the maximal resolving power increases 
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Fig. 2 The maximal mass resolving powerRmaz achieved at 


different number of revolutions and different beam size. The 
lines are to guide the eye. 
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Fig. 3 Potential distribution at the optical axis z in one of 
the mirrors of the MRTOF mass analyzer for the average ion 
kinetic energy of 1500 eV. The position z = 0 corresponds to 
the middle point between the ion mirrors. E; and Ez are two 
regions with approximately constant electric fields. The solid 
and dash lines correspond to the potential distributions when 
the ions are injected-into/ejected-from the analyzer and are 
trapped inside, respectively. 


steadily and then reaches its maximum around N ~ 
150. That means, we can get a very good separation of 
ions in only 3 ms for Artt. 


For every optimization for a specific number of rev- 
olutions we obtained one potential set with the highest 
mass resolving power. Table 1 lists the potentials opti- 
mized for the number of revolutions of N=350. Fig. 3 
shows the potential distribution along the optical axis 
in one of the mirrors of the analyzer when the ions are 
injected-into/ejected-from the analyzer and are trapped 
inside. It has two regions with approximately constant 


Wen-Xue Huang et al. 


Table 1 Optimized potentials on the electrodes for the ion species of 4° Artt. M1 — M4: mirror electrodes 1 — 4 from outside 


to inside, L: lens electrode, D: drift-tube. 


Electrode M1 M2 


M3 M4 L D 


Potential (V) 1373.4 867.9 


962.1 


—1046.2 —2436.7 459.8 


electric fields [25,26], EB, < Ey. The ions turn around 
in the weaker constant electric field F}. 

Applying the potentials in Table 1 to the electrodes, 
we calculated the temporal width at the detector plane 
with respect to the number of revolutions. Fig. 4 shows 
the calculation results of the total TOF, the temporal 
width ATOF, and the mass resolving power R as func- 
tion of number of revolutions. Increasing the number of 
revolutions, the total time-of-flight almost increases lin- 
early, and the initially high temporal spread decreases 
until a minimum, only 61 ns at the number of revo- 
lutions of N ~ 100, and then increases. The maximal 
resolving power R = 3.2 x 104 are achieved at 350 laps, 
corresponding to a TOF of 7.0 ms. 


4.2 Effects of potential inaccuracy 


Fig. 5 shows the relative variation of TOF and mass re- 
solving power as functions of relative variation of bias 
potentials as determined by our optimization code. For 
the total TOF, the behavior from the mirror electrode 
M3 acts differently from other electrodes, and the largest 
variation is from the mirror electrode M1. The maximal 
achievable mass resolving power decreases as the poten- 
tials applied on the electrodes vary from the ideal ones 
as expected. Except M1 and M2, the potential varia- 
tions from all other electrodes (M3, M4, L and D) al- 
most have no effect on the final resolving power. Those 
from the mirror electrode M1 has the biggest effect, 
the resolving power decreases 22% from 3.2 x 104 to 
2.5 x 104 as the relative variation of the potential by 
500 ppm. Thus we have to assure the inaccuracy of the 
potentials applied to the mirror electrodes M1 and M2 
to be less than 50 ppm or preferably 20 ppm. It is very 
easy to understand that the potentials on the electrodes 
M1 and M2 form the mirrors and the ions turn around 
at the position between M1 and M2. 


4.3 Effects of different beam size 


Fig. 2 also shows the maximal mass resolving power 
achieved at different beam sizes. For the same num- 
ber of resolutions, we keep all the parameters of the 
initial beam conditions except the beam size. For all 
three beam sizes, the maximal resolving powers increase 


steadily and then reach their maxima as increasing the 
number of revolutions, but the reaching point delays 
from N ~ 150 to ~ 200 and ~ 300 for the beam size 
reducing from 2.8 mm to $1.4 mm and 0.3 mm, re- 
spectively. 

For different beam size, the maximal resolving power 
increases as the beam size reduces. It reaches 3.2 x 104 
for beam size of 2.8 mm, increases to 4.4 x 104 for 
1.4 mm, and to 5.6 x 104 for 60.3 mm. The smaller 
beam size benefits the resolving power, thus we should 
inject ion beam with a smaller beam size into the mass 
analyzer to achieve a higher mass resolving power. 


4.4 Comparison between the mirror-switching and the 
in-trap-lift modes 


Since we have obtained the optimal parameters for the 
same mass analyzer operating both in the mirror-switching 
and the in-trap-lift modes, we can compare them di- 
rectly. In our optimization, we used the same initial 
beam conditions described above, the initial beam size 
is 2.8 mm diameter and the initial average kinetic en- 
ergy of the ions is 1500 eV. The outcome differences 
only come from the different operating mode. 


In the mirror-switching mode, the maximal resolv- 
ing power has been achieved to be 1.3 x 10° with a total 
time-of-flight of 6.5 ms for the ion species of Artt; 
but in the in-trap-lift mode, it is only 3.2 x 10* with 
a total time-of-flight of 7.0 ms. It seems that the mirror- 
switching mode is much better than the in-trap-lift mode 
if we ignore the financial costs of the power supplies 
and the switches to switch the potentials applied on 
the different electrodes. The difference in mass resolving 
power may come from the difference of the kinetic en- 
ergy of the ions travelling in the trap, which is 1500 eV 
in the mirror-switching mode and 1040 eV in the in- 
trap-lift mode in our case. We may promote the mass re- 
solving power in the in-trap-lift mode by increasing the 
kinetic energy of the ions in the trap. A paper about the 
relationship between maximal resolving power, kinetic 
energy in the trap, mass of ion species, beam size, and 
so on will be prepared after more very time-consuming 
calculations. 
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Fig. 4 Calculated results of the TOF (a), the temporal width ATOF (b), and the mass resolving power R (c) as function of 
number of revolutions. 
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Fig. 5 Calculated relative variation of TOF (left) and mass resolving power R (right) as functions of relative variation of bias 
potentials as determined by optimization code. M1 — M4: mirror electrodes 1 — 4 from outside to inside, L: lens electrode, D: 
drift-tube. 
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5 Summary 


Our previous method to design a multiple-reflection 
time-of-flight mass analyzer has been updated. Now 
constraints can be applied on the potentials. The method 
can be used to search the optimal potentials applied on 
the electrodes in the MRTOF mass analyzer operating 
not only in the mirror-switching mode but also in the 
in-trap-lift mode. 

By using this updated method, the potentials ap- 
plied on the electrodes of our MRTOF mass analyzer 
operating in the in-trap-lift mode have been searched 
and optimized. With a beam size of 2.8 mm diameter 
and an initial average ion kinetic energy of 1500 eV, the 
maximal mass resolving power has been achieved to be 
3.2 x 104 with a total time-of-flight of 7.0 ms for an ion 
species of 40 Artt., The maximal resolving power can 
reach up to 5.6x104 for a beam size of 0.3 mm diame- 
ter. The simulation also shows that the variations of the 
potentials applied to the mirror electrodes M3-M4, the 
lens electrode L and the drift-tube D almost have no 
effect on the final resolving power. Crucially, the inac- 
curacy of the potentials on the mirror electrodes M1-M2 
needs to be less than 50 ppm or preferably 20 ppm. 
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